Masoud Taefi Feijani; Saeed Azadnejad
Abstract
Extended Abstract
1- Introduction
The present study primarily sought to present a new FCD model to eliminate two limitations of the initial FCD model.These limitations included the fact thatimplementing the initial FCD model for sensors without a thermal bandwas not possible, sincethe model ...
Read More
Extended Abstract
1- Introduction
The present study primarily sought to present a new FCD model to eliminate two limitations of the initial FCD model.These limitations included the fact thatimplementing the initial FCD model for sensors without a thermal bandwas not possible, sincethe model took advantage of a combination of shadow index and thermal index to detect black soil and calculate advanced shadow index.To overcome this limitation, we replaced thermal index with NLI and GNDVI indices, and combined shadow, BI, NLI and GNDVI indices to detect black soil and calculate advanced shadow index.Defining a global threshold for thermalindex used fordetectingblack soil was the next limitation of initial FCD model.Due to variations in regional climate and temperature, selecting a global threshold for the whole scene does not seem logical.Thus, a local thresholding process was used to define a threshold level for BI, NLI and GNDVI indices.In this regard, the study sitewas divided into 14 sections and an appropriate threshold wasselected for each section.A digital elevation model was also used to define a specific threshold level for forests in flat areasand elevated areas.
2- Materials & Methods
2.1 Study area and dataset description
The present study was performed within the basin of the Caspian Sea.Drainage basin is considered to be a standard unit ofstudy in environmental studies and thus due to the applied nature of the present study, the Caspian Basin was selected as our study site. In this study, a new FCD model was implemented for data collected from Landsat 5(1366) and Landsat 8 (1396).
2.2 Proposed approach
In the present study, an improved FCD model was obtained by adding two steps to the initial FCD model. In the following paragraph, these two steps will be explained.
2.2.1 Removing thermal index
The first limitation of the initial FCD model lies in the fact that implementing this model for data collected bysensors without thermal band is impossible, because advanced shadow index in the initial FCD model is calculated by combining shadow and thermal indices. Thermal index is only used to separate the shadow of vegetation cover from black soil.In order to overcome this limitationin the improved FCD model, thermalindex is replaced with NLI and GNDVI indices. In this way, black soil and vegetation shadows are separatedusingacombinationofshadow, BI, NLI and GNDVI indices.
The NLI index can be calculated using(1):
(1)
The GNDVI index is also calculatedusing(2):
(2)
2.2.2 Local thresholding
In the initial FCD model,black soil identification and shadow index improvement (advanced shadow index calculation) wereperformedusingthresholdingand based on the combination of shadow and thermal indices.In this model, a number is selected as the threshold of the heat index, and shadow index pixels with values less than this threshold are considered as black soil.Obviously, it is practically impossible to define a threshold and calculate advanced shadow index for large scale areas.
Localthresholding is a much more accurate method of thresholding, which is also used in the improved FCD model.In this method, image received from the study site was divided into 14 sections and a suitable threshold value was selected for BI, NLI and GNDVI indices in each section to calculate advanced shadow index.
Moreover, different thresholds were selected forforests in flat areasand elevated areas.In this regard, digital elevation model of the region was used to separate low-altitude and high-altitude areas.
3. Discussion& Conclusion
Results indicated that the proposed improved FCD model has provided a more accurate estimate of forest canopy density as compared to the initial FCD model.
According to the results, the overall accuracy and kappa coefficient of the initial FCD model were 86.24% and 68.43%, respectively.However, the improved FCD model had an overall accuracy of 96.98% and a kappa coefficient of 92.31% which confirms improved performance of the model.
Moreover, the statistical analysis of changes in the canopy densityindicated that the total area of Hyrcanian forests increased by about 161,963 hectares from 1366 to 1396. This includes an increase ofabout 79, 50 and 33 hectares in Mazandaran, Gilan and Golestan provinces, respectively.
Saeed Azadnejad; Yasser Maghsoudi
Abstract
Extended Abstract
Introduction
Persistent Scatterer Interferometry (PSI) is a technique for detection and analysis of a network of coherent pixels referred to as the Permanent/Persistent scatterer (PS) which have high phase strength over long time periods. This technique has been widely used ...
Read More
Extended Abstract
Introduction
Persistent Scatterer Interferometry (PSI) is a technique for detection and analysis of a network of coherent pixels referred to as the Permanent/Persistent scatterer (PS) which have high phase strength over long time periods. This technique has been widely used by the scientific community to measure the displacement related to thesubsidence/uplift, landslide, tectonic, and volcanoes. As the density and quality of PS pixels are important factors in PSI algorithms, the concept of polarimetric optimization in the PSI algorithms was proposed to improve the number of PS pixels. The recent launch of radar sensors operating with a polarimetric configuration can help improvingthePS-InSAR analysisby increasing the PS density. Therefore, the combination of thepolarimetric and interferometric techniques helpsimprove the PSI techniques, especially in non-urban areas which suffer from lack of the PS density. In this study, we investigated how the contribution of the S1A and TSX data in the PSI analysis could lead to the improvement of the results of the PSInSAR algorithm. Indeed, the main objective of this paper is to illustrate the capability of each dataset for improving the polarimetric optimization results.
1. Materials & Methods
2.1
The proposed method was tested using a dataset of 40 dual-pol SAR data (VV/VH) acquired by Sentinel1-A between February 2017 and May 2018 and 20 dual-pol SAR data (HH/VV) acquired by TerraSAR-X betweenJuly 2013 and April 2014.
2.2 Polarimetric SAR Interferometry
The general principle of polarimetric SAR interferometry was proposed by (Cloude & Papathanassiou, 1997) for the first time. The scattering matrix S represents the polarimetric information associated with each pixel of the image. Considering the monostatic configuration, the scattering matrix S is defined as follows:
(1)
Where and are co-polar channels, is the cross polar channel. This matrix can be represented with the target scattering vector as:
(2)
Where, is the transposed operator. The Pauli vector for the dual-pol data (HH/VV) of the TerraSAR-X sensor, is written as :
(3)
Similarly,the Pauli vectorfor the dual-pol data (VV/VH) of theSentinel1-A sensorcan be expressed as:
(4)
In order to generate scattering coefficient μ, projecting the scatteringvector on the projection vectorwould be sufficient:
(5)
Where is thelinear combination of the elements of matrix S, i is the correspondent of the 2 images, and * represents the conjugate operator. The projection vectorfor the dual-pol data isdefined as:
(6)
Where, and are two real parameters whose ranges are finite and known and are related to the geometrical and electromagnetic properties of the targets. In our research, the main purpose of the polarimetric optimization is to find theoptimum projection vector, in a 2-dimensional search space, and
2.3 Amplitude Dispersion Index Optimization
Substituting (5) into (7), the ADIfor the polarimetric case () can be expressed as follows:
(7)
(8)
According to (6), the polarimetric optimization problem isreduced to finding a suitable and in a finite and known range,so that (8) is minimized.
2. Results & Discussion
The results showed that the proposed method improved the performance of the PSInSAR algorithm in two terms of phase quality and density of the PS pixels. Compared with the VV channel, , the number of PSC and PS pixels increased about 2 and 1.7 times In S1A data, using the ESPO method while, compared with the normal channels like HH and VV, the number of PSC and PS pixels in ESPO method increased about 3.5 and 3 times in TSX data.Based on these results, the optimization methods are more effective in improving the quality of the PSC densitythan in increasing the number of PS pixels. This is mainly because the employed optimization is based on minimizing ADI criterion which is used in the PSC selection. Moreover, ESPO method has been more successful for TSX data compared to the S1A data. This result is due to the higher capability of the TSX data in creating more diverse scattering mechanisms and hence identifying more optimum scattering mechanism compared to S1A data. We also investigated the effect of polarimetric optimization in increasing the PS density in urban and non-urban areas. The experimental results showed that the method succeeded to significantly increase the final set of PS pixels in both urban and nonurban areas.
3. Conclusion
The results show that the optimization methods have been more successful in the improvement of PS density for the TSX data compared to the S1A data. This result is due to the higher capability of the TSX data in creating more diverse scattering mechanisms compared to the S1A data. In summary, thanks to the polarimetric data, it is possible to exploit a larger number of pixels compared with the single polarization case.